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A theory for thermomechanical behavior of homogeneous DNA at thermal equiUbrium predicts 
critical temperatures for denaturation under torque and stretch, phase diagrams for stable B-DNA, 
supercoiling, optimally stable torque, and the overstretching transition as force-induced DNA melt- 
ing. Agreement with available single molecule manipulation experiments is excellent. 
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DNA is a highly refined nanomechanical object. The 
interplay between strong covalent bonds of the backbone 
and weak hydrogen interactions between bases the 
thermal bath in which it is immersed, and the proximity 
of physiological conditions to the denaturation temper- 
ature, make DNA highly and non-linearly responsive to 
mechanical and thermal changes, and render any solely 
mechanical approach inviable. This is critical for nan- 
otechnology, where DNA is the basis for novel materi- 
als [2j , but understanding double helix thermomechanics 
would also illuminate biology, where enzymes involved in 
replication and repair are viewed as molecular motors. 

In the past twenty years, direct single molecule manip- 
ulation [21 13] has revolutionized our understanding of key 
aspects of DNA, revealing new couplings and transitions 
between different structures, whose nature and forms, 
however, are still speculative. When a few micrometer 
long strand of DNA is stretched to a tension of the order 
of pico-Newtons (pNs) to avoid formation of plectonemes, 
sharp transitions are activated at positive and negative 
torques, while an overstretching transition is observed 
for DNA under tension of ~ 60 pN at zero torque [5]. 
Tentative tension-torque phase diagrams for the stabil- 
ity of B-DNA, and various phenomenological theories, 
have been proposed to explain these effects. The ro- 
bustness of the Peyrard-Bishop-Dauxois approach (PBD) 
[BHS] has been corroborated by Cocco and Barbi [3 ES] : 
they incorporated torque and successfully reproduced de- 
naturation by unwinding. However, these recent studies 
do not include tension, do not explain denaturation at 
overwinding, and do not provide phase diagrams in the 
tension-torque plane. Also, although much simpler than 
the molecular structure they describe, their complexity 
still calls for numerical treatments, and cannot offer an- 
alytic equations to more easily guide experiments. 

We address these issues by modeling the steric depen- 
dence of the base bond and the effect of tension and 
torque in a way suitable for elimination of the angular de- 
grees of freedom via integration of the partition function, 
and obtain an effective energy for a PBD model which in- 
corporates temperature and external loads. We compute 
the phase diagram for B-DNA and the dependence of 
supercoiling on torque, tension and temperature at crit- 
icality. Finally, we propose simple algebraic formulae for 



the observables. 

In our model i labels nucleotides separated by a dis- 
tance a along the DNA backbone (Fig. 1), Xi is the length 
of the z*'' base's bond, oji = {9i+i — 6i-i)/2 — His the an- 
gular shift between nucleotides along the backbone, and 
6i is their angular coordinate. As in the Cocco-Barbi 
models the two strands of DNA are assumed symmetri- 
cal even when open, il denotes the natural pitch of the 
helix, and uji describes deviations from equilibrium. The 
potential energy of the system is E = a Ei with 

E,^'l^+'^{u;, + nf + [x{x,) - 1] Vioj,) (1) 

the sum of: a stacking potential (Axi = Xi+i— Xi), in har- 
monic approximation for simplicity (see discussion later) ; 
an elastic term which restores the = Qi angular con- 
figuration for open strands; and a square well potential 
for the hydrogen bond between bases [x(a;) is a step func- 
tion which is for Q < x < Xc and 1 for x^ < x, where 
Xc is a length associated with the hydrogen bond] , whose 
depth V depends on the angle w^. Because of a com- 
plex combination of hydrophobic, tt-tt, and dipolar in- 
teractions, the bonding of opposite bases is responsible 
for DNA's pitch. Therefore V^(cli) is not symmetric but 
rather V^'(O) = since DNA is in equilibrium at — 0. 
Also, [I = —V"{Q) must be positive: indeed, i/ + ^ being 
the torsional rigidity of the joined double helix and i/ of 
the (much softer) open strands, we have fj. ^ ly > (we 
shall see that /i/i^ ^ 10^). We can now expand V as 

V{uj) ~Vo + vVLuj - ^/iw^. (2) 

Our choice of a square well potential and separation of 
variables x^, in ([T]) keeps the model analytically solv- 
able: no substantial changes in the thermodynamics arise 
from surrendering it in favor of, e.g., the Morse poten- 
tial 8J. We show elsewhere how to introduce a smooth 
potential, desirable for dynamics and other studies [TT]. 

Now consider the external loads. The torque F is in- 
corporated in ([l| via a term TUJ^a (for dimensional 
convenience F = to). Tension is more subtle: the to- 
tal stretch hi = Y,i x{xi)h° + J^i ixi^i) - 1] from 
both closed {hf) and open {h°) DNA sections can only 
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FIG. 1: Schematics of our DNA model showing degrees of 
freedom and geometrical parameters (Left), closed double he- 
lix (Middle), and open strands twisting with bases pointing 
outwards (Right). 
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arise from winding/unwinding, since the backbone is ef- 
fectively inextensible. While elongation due to a change 
in pitch for the joined double helix is trivially (h'^+ho)'^ — 
hi - R^iujf + 2a;, n) {R is the radius of the DNA helix, 
hf) < a the vertical distance between nucleotides), for 
open sections assumptions are necessary. We speculate 
that the two backbones twist with bases pointing out- 
wards as in P-DNA (Fig. 1), and therefore the length of 
openings responds to winding: = — r'^{uJi + il)^, 
where r < i? is an effective diameter for the backbone. 
Below, we will expand these expressions to second or- 
der around their stable configurations (wi — for closed, 
LOi = — r2 for open sections). All the quantities k, 
H, Vb, T, and / have the dimension of a force. 

Equilibrium thermodynamics is implemented by inte- 
gration of exp —f3{E — ^ ■ TUJiQ — fhi) over {xi} and 
{uJi}. Since everything is quadratic in Wi, we can ex- 
press the product of the gaussian integrals in the angular 
variables and obtain the partition function 



Z = 



-/3a 



-+V{xi} 



(3) 



for an equivalent PBD model, whose effective potential 



V{xi) = [x{x,) - 1] 
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incorporates explicitly the effect of the external torque 
r — AT, and also of /, through the tension-increased 
torsional rigidities fl = fi + mf, v ^ v + nf, and the 
pitch under tension il ~ il — i^^f- Also, the depth of 
the effective potential in the absence of external torque, 

weakened by tension / and entropically by temperature 
T. The term A = -fln^ - - vf - nr - :^T^ 

from while irrelevant for the phase diagram, must 
be kept when computing supercoiling. L — Na is the 
length of the DNA. There are then purely geometrical 
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10 A) is the radius of the DNA 



FIG. 2: Predicted critical surface for denaturation of B-DNA 
in (T), torque (F) and tension (/) obtained from ([7|. Solid 
lines on the surface correspond to critical lines at fixed torque 
r plotted in Fig. [5] 



molecule, ho (~ 3.4 A) the elevation between consecutive 
nucleotides, a (~ 7 A) their distance along the backbone, 
and = 27r/10 the pitch of DNA: these are established 
geometrical values for B-DNA, but our formalism works, 
mutatis mutandis, for A- and Z- forms. 

The expression for Z in ([s]) is exact within our model, 
and the resulting PB problem (or PBD if anharmonic- 
ity is included) is amenable to numerical transfer matrix 
treatment. Here we will proceed analytically. In the con- 
tinuum limit and neglecting an irrelevant equipartition 
factor, Z in ([3]) can be written as 



Z oc e-^^^ Tr e"^^, 

proportional to the trace of the operator 
1 



H 
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d', + pvix,). 



(5) 
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Torque, tension and temperature enter the potential (|4|, 
and the critical surface corresponds [B] to the disappear- 
ance of the bound state for H, or 



Ht — T — wf H 

2(i> + /i)i> ~ ~ 



1 

2a 



A* 



TdI-TI] =0, 



(7) 



where e — ^kx^jii'^ lias the dimension of an energy, / — 
ln[(z/ + /i)/i^] and I = ln[(i> + /i)/i>] are pure numbers, and 
r|, = eXcSVo - i^ln^ - is the denaturation 

temperature in the absence of torque or tension. 

Figure [2] shows our predicted critical surface of stable 
B-DNA in the space of F, / and T . There are only a few 
parameters to fit: an accepted value for the denaturation 
temperature is To = 350 K [T^; we show below from 
data on torque-winding experiments that [i = 10'^ pN. 
Choosing the remaining three parameters as tx^ = 45 
pNA, = 24 pN and n = 0.3, provides a remarkably 
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FIG. 3: Predicted critical lines for denaturation at different 
temperatures. The region enclosed by each line corresponds 
to stable B-DNA at that temperature. Points are experimen- 
tal data from single molecule manipulation [5], square corre- 
sponds to the well known overstretching transition (F = 0, 
/ = 60 pN). Td = 350 K = 77 °C, the denaturation tempera- 
ture, is critical at zero external load: but even at temperatures 
of Td or higher (inset), B-DNA can be stable in an interval of 
applied positive torque. Solid (dotted) straight line indicates 
Tm ~ avQ. (its approximation Tm ~ 11 -1-0.092/), the middle 
point between critical torques, which also corresponds to the 
highest critical temperature at given tension. 

good fit for 10 experimental data points (20 numbers) [S] 
in the / vs. T phase diagram of Fig. [3] The skewness 
of the critical lines can be quantified: from ([t]) one finds 
for r„j — (r+ +r^)/2, the middle point between critical 
torques at given tension, 

Tm = a^—^ — D Cl c:i ai> n + f ( nfl ), (8) 

independent of temperature. As ([s]) reveals, skewness re- 
sults from backbones twisting under torque: it would be 
erroneously negative for r = 0. With the above parame- 
ters, 

r„/ (pN nm) = 11 + 92 X IQ-^ // (pN) . (9) 

Equations ([Sj |9| along with experimental data for DNA 
that place melting at / = 15 pN, for r+ = 34 pN nm, 
T- = -10 pN nm and at / = 60 pN, for r+ = 33 pN 
nm [S] , predict melting at zero torque and tension / = 60 
pN (blue square in Fig. [3| in good agreement with the 
observed overstretching transition. Our analysis implies 
a force-induced melting [13 rather than a transition to 
a double helix with distortions [14 . A positive torque 
stabilizes DNA even at temperatures above denaturation 
(Fig[3] inset), a property exploited by thermophile bacte- 
ria living at high temperatures [1]. As expected, negative 



FIG. 4: Predicted critical temperature as function of torque 
for stretched DNA at different tensions (0,15, 30, 45, 60, 
75, 80, 85 90 pN). Points are experimental data from sin- 
gle molecule manipulation fE\ performed at 23 °C. Squares 
correspond to the overstretching transition (F = 0, / = 60 
pN) at 23°C and denaturation (F = 0, / = 60 pN) at 
Tn = 77 °C. The maximal critical temperature for a given 
stretch is achieved under external torque Fm(/) given by ([8|, 
and is 91 °C at zero tension and Tm = 11 pN nm. 

torque destabilizes DNA, a mechanism exploited in biol- 
ogy for DNA opening and replication. 

Under a given stretch, the critical temperature in- 
creases (decreases) under positive (negative) torque, as 
shown in Fig. [4] From ([7| we see that T'^ is maximized 
at r = r,„(/), which therefore induces the most sta- 
ble configuration at any temperature, for a certain ten- 
sion /. The highest temperature B~DNA can with- 
stand is achieved under zero tension and positive torque 
Tm = avil = 11 pN nm, and corresponds to — 91 °C. 

DNA supercoils for packaging inside cells. This corre- 
sponds to small positive or negative torques; in that bio- 
logical regime, the critical temperature decreases mono- 
tonically with an applied tension, in fact linearly, with 
slope independent of the applied torque or degree of su- 
percoiling. Yet, for torques larger than about 32 pN nm, 
a regime accessible to single molecule manipulation ex- 
periments and potentially useful in nanotechnology, the 
maximal critical temperature corresponds to a non-zero 
tension, suggesting that at low temperatures and large 
torques DNA can be stabilized by tension (Fig.[5|. 

The average change in pitch is given by (w) = 
i^^/3^^9r In Z = T/v — Vl — TdreB (es is the bound state 
of the hamiltonian ([6]), whereas the remainder comes 
from A). At the critical point, we find 

\l ai' \ it / u + v au 
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FIG. 5: Predicted critical temperature as a function of tension 
for DNA twisted under different torque. Numbers on curves 
denote the external torque in pN nm. The maximum critical 
temperature at any given tension (dotted curve) is reached at 
external torque rm(/) from (|8|. While for negative torque 
the curves show a monotonic and linear decrease (inset), for 
positive torque larger than 32 pN nm (dashed curve) curves 
become non-monotonic: for low temperatures and large ap- 
plied torque, B-DNA is stabilized by tension. 



where ujd is the small negative unwinding (typically 
ud — —10^'' rad) at denaturation (T = Tq, F = 0, 
0). Equation (10) can be used to fit, from experi- 



/ _ 
mental data [5], ^ = 10^ pN. When neglecting small non- 



linear corrections, ( 10 1 is also a good approximation for 



the torque vs. pitch curve away from criticality. While 



the derivation of ( 10 1 will be presented elsewhere fJJ], we 
note here that a well potential on an infinite half line in 
^ would generate no contribution from the bound state 
to (w) at criticality, and thus give (wc) = ^ — fl, the 
expected behavior for separated strands. The disconti- 
nuity at criticality in (w) seen in experiments therefore 
comes from the finiteness of lateral entropy: when bases 
open in a bubble, they cannot move infinitely far from 
each other: as tension and torque are applied and back- 
bones twist, the radius of DNA provides a natural barrier 
for large x. This suggests that in biology, where, unlike in 
manipulations under tension, larger distances are avail- 
able between bases in larger openings, the average change 
in pitch (w) might access larger values before complete 
strand separation occurs. 



Equation (10 1 does not predict the tiny anomalous 



overwinding under tension [TS] . We believe this is a con- 
sequence of neglecting the bending modes, whose ther- 
mal fluctuation weakens the base's bonds. As tension 
straightens DNA, it might reduce this effect by collaps- 



ing openings and thus overwinding the structure. Con- 
versely, torque-induced overwinding reduces the fraction 
of open bases, thereby stiffening the bending modes and 
reducing the amplitude of their thermal fluctuations, re- 
sulting in elongation. We will show elsewhere [11] how 
to incorporate these effects into our formalism to predict 
elongation under external loads, and also other phenom- 
ena relevant to biology, such as bending induced by a 
genetic defect, and effects of non-uniform gene sequence. 
Finally, a transfer integral numerical study [TT] of the 
anharmonic version of ([s]) is needed to augment our pre- 
dictions at temperatures closer to denaturation, where 
the effect of nonlinearity in the stacking potential is im- 
portant for the transition order and precursors |16) . 

In summary, a range of predictions and phenomena 
for DNA thermomechanics (critical lines, phase dia- 
grams, supercoiling under loads, optimally stable torque, 
tension-induced stability at high torque) that were inac- 
cessible to previous models, or could only be covered par- 
tially and numerically, have been explained here within 
a unifying framework amenable to analytical treatment 
and further extensions and applications. We are grateful 
to B. Alexandrov and K. Rasmussen (LANL), S. Ares and 
J. Bois (MPIPKS Dresden), C. Matek and A. Louis (Ox- 
ford) for stimulating conversations, to Amity Law (Har- 
vard), Cynthia Reichhardt (LANL), and Paul Lammert 
(PSU) for helping with the manuscript. 
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